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Abstract. Present and future observations of supernova relic neutrinos (SRNs), i.e., 
a cosmological neutrino background from past core-collapse supernova explosions, 
potentially give us useful information concerning various fields of astrophysics, 
cosmology and particle physics. We review recent progress of theoretical and 
observational studies of SRNs, particularly focusing on the detectability and also on 
implications for cosmic star formation history and neutrino physics. 
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1. Introduction 

A core-collapse supernova explosion is one of the most spectacular events in astrophysics, 
and it attracts a great deal of attention from many physicists and astronomers. It 
also produces a number of neutrinos and 99% of its gravitational binding energy is 
transformed to neutrinos; detection of the galactic supernova neutrino burst by ground- 
based large water Cerenkov detectors, such as Super-Kamiokande (SK) and Sudbury 
Neutrino Observatory (SNO), would provide valuable information on the nature of 
neutrinos as well as supernova physics. In addition, because supernova explosions have 
occurred very commonly in both the past and present universe, tracing the cosmic star 
formation rate (SFR), they should have emitted a great number of neutrinos, which now 
make a diffuse background, i.e., supernova relic neutrinos (SRNs). Involved physics 
in SRN ranges quite widely — from cosmic SFR and supernova physics to neutrino 
properties as elementary particles. Therefore, detecting SRNs or even setting limits 
on their flux can give us quite useful and unique implications for various fields of 
astrophysics, cosmology and particle physics. Detectability of SRNs in various detectors 
and its implications have been discussed in many theoretical papers [T]-|19| from various 
points of view, about which we discuss in detail in the following part of the present paper. 

Flux estimation requires models of neutrino spectrum emitted from each supernova 
explosion and cosmic SFR. Furthermore, it is experimentally established that neutrinos 
mix among different flavours, altering the neutrino spectrum at detectors from the 
original one after neutrinos propagate inside a supernova envelope; this effect should 
now be taken into account appropriately. In addition to a flux estimation as precise as 
possible, a detailed discussion of background events, which hinder the SRN detection, 
is essential; it has been thoroughly studied by Ando, Sato and Totani ^2] and we 
follow their discussion later in this paper. A stringent observational upper limit on 
the SRN flux is obtained by the SK group (2Hj and its value is only a factor of 3- 
6 larger than several theoretical predictions |H] ^OJ , being consistent with (and 
more stringent than) the theoretical upper limit given by Kaplinghat, Steigman and 
Walker jTTJ using conservative models. In order to make the SRN detection more likely, 
a promising method was proposed by Beacom and Vagins [21 J. Their basic idea is to 
dissolve gadolinium trichloride (GdCbj) into the water Cerenkov detectors, which greatly 
reduces the background events if it is applied to the currently working or proposed future 
detectors such as SK, Hyper-Kamiokande (HK) and Underground Nucleon Decay and 
Neutrino Observatory (UNO). Therefore, we are now at an exciting stage, where SRNs 
would soon be detected actually, and be used to obtain several implications for various 
fields of astrophysics as a unique and complementary method to usual observations of 
the light. 

Since we have a promising prospects of the SRN detection, it is timely to discuss 
the SRN potential for probing the universe or particle physics. Cosmic SFR is one 
of such possibilities. The most popular method inferring SFR is to use the galaxy 
luminosity function of rest-frame ultraviolet (UV) radiation 22 3.0) • m addition to far- 
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infrared (FIR) / sub- millimeter dust emission jSTJ E2] and near-infrared (NIR) Ha line 
emission [SSI [SHI - in these traditional approaches, however, there are a fair number of 
ambiguities when the actual observables are converted into the cosmic SFR [HI]. The 
most serious problem is that the effects of dust extinction are non-negligible, especially 
for rest-frame UV observations. For the approach using SRNs, on the other hand, we are 
not troubled by such a problem, because neutrinos are completely free of dust extinction. 
This point is the same with observations in the sub-millimeter wave band; however, 
neutrinos are emitted directly from stars, whereas sub-millimeter radiation comes from 
dust and is an indirect process. Another advantage is that supernovae are directly 
connected with the death of massive stars with M > 8M Q , whose lifetime is expected to 
be very short compared with the Hubble time-scale Hq 1 . It enables direct inference of 
the cosmic SFR assuming the initial mass function (IMF), not bothered with uncertainty 
concerning the galaxy luminosity fuction, which sometimes causes difficulty in the case 
of rest-frame UV observation [37] . Motivated by all of these reasons, many researchers 
have investigated the dependence of the SRN flux on various models of cosmic SFR or 
to what extent the SFR can be probed from the future SRN observations. Totani, Sato 
and Yoshii [H] gave flux estimation using several theoretical models of galaxy evolution. 
After the bursting release of the observational SFR data points since the pioneering 
study by Madau et al 23., the SRN flux has been estimated using these SFR results 
[01 HH1 E2] • More recently, theoretical SRN flux calculations have been compared with 
the observational upper limit by SK [2U] . Fukugita and Kawasaki [13 j used the limit 
to obtain constraints on the cosmic SFR; Strigari et al [IH] adopted the latest SFR 
model based on the Sloan Digital Sky Survey 38 and concluded that their median SRN 
flux is slightly below the current SK upper limit. Following the proposal to dissolve Gd 
into detectors [2T], Ando [T7j investigated potential performance of the Gd- loaded water 
Cerenkov detectors such as SK, HK and UNO, using the Monte Carlo (MC) simulation, 
especially focusing on how reliably an assumed SFR model can be reproduced from the 
SRN observation. 

Neutrino properties can also be probed by the SRN observation in principle. It 
is pointed out that the resulting SRN signal depends on neutrino oscillation models 
[Tl] : particularly if the mass hierarchy is inverted, i.e., the first mass eigenstate that 
most strongly couples to electron flavour is heavier than the third state, m\ > 7713, and 
the value of #13, which has not been well constrained yet, is sufficiently large, then the 
flavour conversion during propagation inside the supernova envelope could considerably 
change the neutrino signal. Another interesting possibility is that a stringent constraint 
can be set on neutrino decay models from the SRN observation as first pointed out by 
Ando [12] and succesively studied by Fogli et al [18]. Non-radiative neutrino decay can 
be induced by the interaction between neutrinos and massless or very light particles 
such as Majoron. The strongest lower limit to the neutrino lifetime-to-mass ratio is 
obtained from the solar neutrino observation [HH] [H] and meson decay [12] 01] to be 
r/m > s/eV. SRN observation, on the other hand, is sensitive to ~ 10 10 s/eV, 
which is many orders of magnitude larger than the current limit. This is because of 
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the much longer baseline between the Earth and cosmological supernovae, from which 
SRNs are emitted, than 1 AU in the case of solar neutrino observations. Neutrino decay 
could also change the detected signal from high-energy astrophysical objects |25] or the 
galactic supernova explosions [IE] [H], an d could alter usual discussions on the early 
universe and structure formation jlHj as well as on supernova coolings |HH]-|52j: but we 
stress that in principle we can obtain the most stringent constraints on the neutrino 
lifetime from the SRN observations. 

This paper is organized as follows. We first introduce formulation for calculating the 
SRN flux in sectional Models of cosmic SFR, supernova neutrino spectrum calculated 
by various groups are given in the same section. Neutrino oscillation inside the supernova 
envelope is briefly described there. SRN flux and event rate at ground-based detectors 
calculated with our reference models are given in section El Section |U is devoted to a 
detailed discussion on background events against the SRN detection and detectability at 
various detectors. Current observational limit by SK, which is obtained by a statistical 
argument including background events, is briefly summarized also in the same section. 
We discuss implications of the SRN observation for cosmic SFR and neutrino properties 
such as oscillation and decay in sections El and El respectively. Finally, we conclude 
this paper by giving brief summary in section [7| Throughout this paper, we only 
consider electron anti-neutrinos (P e ) at detectors because this kind of flavour is most 
easily detected at the water Cerenkov detectors, on which we mainly focus. 



2. Formulation and models 



2.1. Formulation 

The present number density of SRNs (u e ), whose energy is in the interval E v ~ E u + dE u , 
emitted in the redshift interval z ~ z + dz, is given by 

dn v {E v ) = Rbs(z)(1 + zffdz ^^- dE' v {\ + z)~ 3 

= R S v(z)^dz^^(l + z)dE„, (1) 

where E' v = (1 + z)E v is the energy of neutrinos at redshift z, which is now observed 
as E u ; Rsn(z) represents the supernova rate per comoving volume at z, and hence the 
factor (1 + z) 3 should be multiplied to obtain the rate per physical volume at that time; 
dN u /dE u is the number spectrum of neutrinos emitted by one supernova explosion; and 
the factor (1 + z)~ 3 comes from the expansion of the universe. The Friedmann equation 
gives the relation between t and z as 
dz 

— = -H (l + z)y/n m {l + zf + Q A , (2) 

and we adopt the standard ACDM cosmology (fi m = 0.3, = 0.7, and H = 
70 h 70 km s" 1 Mpc _1 ).§ We now obtain the differential number flux of SRNs, dF u /dE u , 

§ Although we use the specific cosmological model here, the SRN flux itself is completely independent 
of such cosmological parameters, as long as we use observationally inferred SFR models; see their 



Supernova relic neutrinos 



5 



using the relation dF v /dE v = cdn u /dE v : 



dF u c f z ™ , dN u {E' v ) dz 



where we assume that gravitational collapses began at the redshift z max = 5. 
2.2. Models for cosmic star formation rate 

As our reference model for the SFR, we adopt a model that is based on recent progressive 
results of rest-frame UV, NIR Ha, and FIR/sub-millimeter observations; a simple 
functional form for the SFR per unit comoving volume is given as 



i(\ n^u exp(3.4z) y/Q m (l + z)* + Sl A 3 

V>,(z) = 0.32/Ao 7— r . r , M yr Mpc , 4 

exp(3.8z) +45 (1 + 

where /* is a factor of order unity, which we illustrate below. Figure ^ shows the 
SFR if)*(z) with the various data points from rest-frame UV [22H23E01, Ha line 

and FIR/sub-millimeter E2j observations; these data points are not corrected 
for dust extinction. In the local universe, all studies show that the comoving SFR 
monotonically increases with z out to a redshift of at least 1. Although there is such a 
general observational tendency, even the local (z = 0) SFR density is far from precise 
determination; it ranges fairly widely as V*(0) =(0.5-2.9) xKrt 70 M yr" 1 Mpc~ 3 
|54j . For this reason, we introduce the correction factor /*. Our reference model (|4*]l 
with /* = 1 is consistent with mildly dust-corrected UV data at low redshift; on the 
other hand, it may underestimate the results of the other wave band observations. In 
fact, it predicts the local SFR value of 0.7 x 10~ 2 hy M Q yr -1 Mpc~ 3 , which is close to 
the lower limit of the estimation by Baldry and Glazebrook [54J, but we stress that the 
SRN flux given below using this model can be simply applied to the other cases if we 
adjust the correction factor /*, and would be quite general. Calculations by Strigari et 
al [IB] are based on the SFR model giving local value of 1.6 x 10 -2 M yr -1 Mpc -3 , 
which corresponds to /* = 2.3 in our notation. Although the SFR- z relation generally 
tends to increase from z — to ~ 1, behaviours at the higher redshift region z > 1 are 
not clear at all. Ando et al ^2] also investigated the dependence on the several adopted 
SFR models, which were only different at high-redshift regions (z > 1.5); our reference 
model (jlj) was referred to as the "SF1" model there. They showed that the SRN flux at 
E v > 10 MeV is highly insensitive to the difference among the SFR models (owing to 
the energy redshift, as discussed in section EH} ■ 

We obtain the supernova rate (-Rsn(z)) from the SFR by assuming the Salpeter IMF 
(4>(m) oc m -2 ' 35 ) with a lower cutoff around O.5M , and that all stars with M > 8M & 
explode as core-collapse supernovae, i.e., 

/ Q 25 dm m(p(vi ) 
cancellation between equations J3J and (J3J. 



^sn(^) = *SL , —~M Z ) = 0.0122M"V*(4 (5) 
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Figure 1. Cosmic star formation rate as a function of redshift. Data points are given 
by rest-frame UV (open triangles) E3 E0|, NIR Ha (crosses) |S3 El 113, and 
FIR/sub- millimeter (filled diamonds) 132) observations. The solid curve represents 
our reference model given by equation Q with = 1 . The standard ACDM cosmology 
is adopted (fi m = 0.3, Q A = 0.7, H = 70 km s" 1 Mpc" 1 ) . 



Here we assume that the IMF does not change with time, which may be a good 
approximation provided there are no significant correlations between the IMF and 
the environment in which stars are born; extant evidence seems to argue against 
such correlations over the redshift range of interest (z < 2) The resulting local 

supernova rate evaluated with /* = 1 agrees within errors with the observed value of 
-Rsn(O) = (1.2 ± 0.4) x lCT 4 /^ yr -1 Mpc~ 3 (e.g., [HE] and references therein). In fact, 
the totally time-integrated neutrino spectrum from massive stars (> 30M Q ) could be 
very different from the models that we use (and give in the next subsection), possibly 
because of, e.g., black hole formation. However, the conversion factor appearing in 
equation (jSj) is highly insensitive to the upper limit of the integral in the numerator; for 
instance, if we change the upper limit in the numerator to 25M , the factor becomes 
O.OIOMq 1 , which is only slightly different from the value in equation @. 

2.3. Neutrino spectrum from supernova explosions 

For the neutrino spectrum from each supernova, we adopt three reference models by 
different groups, i.e., simulations by the Lawrence Livermore (LL) group [HZ] and 
Thompson, Burrows and Pinto jHEj (hereafter TBP), and the MC study of spectral 
formation by Keil, Raffelt and Janka [SH] (hereafter KRJ). In this field, however, the 
most serious problem is that the recent sophisticated hydrodynamic simulations have not 
obtained the supernova explosion itself; the shock wave cannot penetrate the entire core. 
Therefore, many points still remain controversial, e.g., the average energy ratio among 
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Table 1. Fitting parameters for supernova neutrino spectrum. 



Model 



Mass E 9a E Ul 

(M ) (MeV) (MeV) /3 Pc 




LL |S7] 

tbp ED 



20 
11 
15 
20 



15.4 
11.4 
11.4 
11.9 
15.4 



21.6 
14.1 
14.1 
14.4 
15.7 



3.8 
3.7 
3.7 
3.6 
4.2 



2.2 
2.2 
2.2 
2.5 



1.8 4.9 x 10 52 5.0 x 10 52 



krj inm 



neutrinos of different flavours, or how the gravitational binding energy is distributed to 
each flavour. All these problems are quite serious for our estimation, since the binding 
energy released as v e changes the normalization of the SRN flux, and the average energy 
affects the SRN spectral shape. Traditionally, neutrino spectrum is assumed to have a 
Fermi-Dirac spectral shape with T Pe ~ 5 MeV and T^ x ~ 8 MeV, where v x represents 
non-electron neutrinos and anti-neutrinos, as adopted in many studies including recent 
ones [7j-[H], [EE EE EH] • This approximation is roughly consistent with the LL model, 
although slight difference exists at both high- and low-energy regions [T2"l l57|. The other 
two models (TBP and KRJ) that are more sophisticated, however, are not consistent 
with such simple treatment at all as we describe below or as already discussed in [T7j . 
Because of all these reasons stated above, we believe that these three models from 
different groups will be complementary. 

The numerical simulation by the LL group 57 is considered to be the most 
appropriate for our estimation, because it is the only model that succeeded in obtaining 
a robust explosion and in calculating the neutrino spectrum during the entire burst 
(~ f5 s). According to their calculation, the average energy difference between z/ e and 
u x was rather large and the complete equipartition of the binding energy was realized 
L Uc = Lp c = L„ x , where L Va represents the released gravitational energy as a-flavour 
neutrinos. The neutrino spectrum obtained by their simulation is well fitted by a simple 
formula, which was originally given by KRJ as 



where E v is the average energy; the values of the fitting parameters for the z/ e and z/ x 
spectrum are summarized in Tabled 

Although the LL group succeeded in obtaining a robust explosion, their result has 
recently been criticized because it lacked many relevant neutrino processes that are now 
recognized as important. Thus, we adopt the recent result of another hydrodynamic 
simulation, the TBP one, which included all the relevant neutrino processes, such as 
neutrino bremsstrahlung and neutrino- nucleon scattering with nucleon recoil. Their 
calculation obtained no explosion, and the neutrino spectrum ends at 0.25 s after core 
bounce. In the strict sense, we cannot use their result as our reference model because the 
fully time-integrated neutrino spectrum is definitely necessary in our estimate. However, 
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we adopt their result in order to confirm the effects of recent sophisticated treatments of 
neutrino processes in the supernova core on the SRN spectrum. The TBP calculations 
include three progenitor mass models, i.e., 11, 15 and 20M Q ; all of these models are 
well fitted by equation (jUJ), and the fitting parameters are summarized in tabled The 
average energy for both z/ e and z/ x is much smaller than that by the LL calculation. 
Although we do not show this in table ^ it was also found that at least for the early 
phase of the core-collapse, the complete equipartition of the gravitational binding energy 
for each flavour was not realized. However, it is quite unknown whether these trends 
hold during the entire burst. In this study, we adopt the average energy given in tabled 
as our reference model, while we assume perfect equipartition between flavours, i.e., 
L Po = L Vx = 5.0 x 10 52 erg. 

In addition, we also use the model by KRJ. Their calculation did not couple with 
the hydrodynamics, but it focused on the spectral formation of neutrinos of each flavour 
using an MC simulation. Therefore, the static model was assumed as a background of 
neutrino radiation, and we use their "accretion phase model II," in which the neutrino 
transfer was solved in the background of a 150 ms postbounce model by way of a 
general relativistic simulation. The fitting parameters for their MC simulation are 
also summarized in table ^ Unlike the previous two calculations, their result clearly 
shows that the average energy of z/ x is very close to that of z/ e . It also indicates that 
the equipartition among each flavour was not realized, but rather L Ue ~ L Dc ~ 2L;, x . 
However also in this case, since the totally time-integrated neutrino flux is unknown from 
such temporary information, we assume perfect equipartition, L Dc = L v ^ = 5.0 x 10 52 
erg, as well as that the average energies are the same as those in table ^ 

2.4- Neutrino spectrum after neutrino oscillation 

The original z/ e spectrum is different from what we observe as z/ e at Earth, owing to 
the effect of neutrino oscillation. Since the specific flavour neutrinos are not mass 
eigenstates, they mix with other flavour neutrinos during their propagation. The 
behaviour of flavour conversion inside the supernova envelope is well understood, because 
the relevant mixing angles and mass square differences are fairly well determined by 
recent solar [HH1 EI], atmospheric |62| . and reactor neutrino experiments jHSj- The 
remaining ambiguities concerning the neutrino oscillation parameters are the value 
of #13, which is only weakly constrained (sin 2 6*13 < 0.1 [03]), and the type of mass 
hierarchy, i.e., normal (m>i <C m 3 ) or inverted (m>i ^> m 3 ). We first discuss the case of 
normal mass hierarchy as our standard model; in this case, the value of #13 is irrelevant. 
The case of inverted mass hierarchy is addressed in section 16.11 In addition, other 
exotic mechanisms, such as resonant spin-flavour conversion (see [65J and references 
therein) and neutrino decay [H3IIB], which possibly change the SRN flux and spectrum, 
might work in reality, and these topics are also discussed later in sections 16.21 and 16.31 
respectively. 

The produced v e at the supernova core are coincident with the lightest mass 
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eigenstate \>\ owing to the large matter potentials. Since this state V\ is the lightest 
also in vacuum, there are no resonance regions in which one mass eigenstate can change 
into another state, and therefore i> e at production arrives at the stellar surface as V\. 
Thus, the u e spectrum observed by the distant detector is 

dE Pe dE Pl dEp 2 dEp 3 

HC/ol|2 g + (1 _ |c , oin g, (7) 

where the quantities with superscript represent those at production, U a i is the mixing 
matrix element between the a-flavour state and i-th mass eigenstate, and observationally 
\U e i\ 2 = 0.7. In other words, 70% of the original z/ e survives; on the other hand, the 
remaining 30% comes from the other component z/ x . Therefore, both the original i> e 
and v x spectra are necessary for the estimation of the SRN flux and spectrum; since the 
original z/ x spectrum is generally harder than that of the original u e , as shown in tabled 
the flavour mixing is expected to harden the detected SRN spectrum. 



3. Flux and event rate of supernova relic neutrinos 

3.1. Flux of supernova relic neutrinos 

The SRN flux can be calculated by equation © with our reference models given in 
section El Figure Efa) shows the SRN flux as a function of neutrino energy for the 
three supernova models, LL, TBP and KRJ. The flux of atmospheric v e , which becomes 
background events for SRN detection, is shown in the same figure pT> The SRN 
flux peaks at < 5 MeV, and around this peak, the TBP model gives the largest SRN flux 
because the average energy of the original v e is considerably smaller than in the other 
two models but the total released energy is assumed to be the same. On the other hand, 
the model gives a smaller contribution at high-energy regions, E v > 10 MeV. In contrast, 
the high-energy tail of the SRN flux with the LL model extends farther than with the 
other models, and it gives flux more than a order of magnitude larger at E v = 60 MeV. 
This is because the high-energy tail was mainly contributed by the harder component of 
the original neutrino spectrum; in the case of the LL calculation, the average energy of 
the harder component v x is significantly larger than that of the other two calculations, 
as shown in table ^ We show the values of the SRN flux integrated over the various 
energy ranges in table EJ The total flux is expected to be 11-16 /* cm -2 s _1 for our 
reference models, although this value is quite sensitive to the shape of the assumed SFR, 
especially at high- 2: . The energy range in which we are more interested is high-energy 
regions such as E v > 19.3 MeV and E u > 11.3 MeV, because as discussed below, the 
background events are less critical and the reaction cross section increases as cx E 2 . In 
such a range, the SRN flux is found to be 1.3-2.3 /* cm" 2 s" 1 [E v > 11.3 MeV) and 0.14- 
0.46 /* cm" 2 s _1 (E v > 19.3 MeV). Thus, the uncertainty about the supernova neutrino 
spectrum and its luminosity gives at least a factor 2-4 ambiguity to the expected SRN 
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Table 2 


. Flux of supernova relic neutrinos. 




Model 


Redshift range 




Flux (/* cm 2 s 


- 1 ) 


Total 


E v > 11.3 MeV 


E v > 19.3 MeV 


LL 


Total 


11.7 


2.3 


0.46 




< z < l a 


4.1 (35.3) 


1.6 (70.9) 


0.39 (85.2) 




1 < z < 2 a 


4.9 (42.0) 


0.6 (26.3) 


0.06 (14.0) 




2 < z < 3 a 


1.8 (15.1) 


0.1 ( 2.5) 


0.0 ( 0.7) 




3 < z < 4 a 


0.6 ( 5.3) 


0.0 ( 0.2) 


0.0 ( 0.0) 




4 < z < 5 a 


0.2 ( 2.1) 


0.0 ( 0.0) 


0.0 ( 0.0) 


TBP 


Total 


16.1 


1.3 


0.14 


KRJ 


Total 


12.7 


2.0 


0.28 



a Contribution from each redshift range to the total (0 < z < 5) value are shown in 
parentheses as percentages. 



flux in the energy region of our interest. 



1 1 I r 
. LL 



TBP 

KRJ 

Atm— v. 




10 20 30 40 50 60 

Neutrino Energy [MeV] 



1 1 I 

. Total 
0<z<l 




10 20 30 40 50 60 
Neutrino Energy [MeV] 



Figure 2. (a) SRN flux in units of /* cm -2 s _1 MeV -1 calculated with three reference 
models of original neutrino spectrum: LL, TBP and KRJ. The flux of atmospheric 
neutrinos [23 EH] is also shown for comparison, (b) The same as (a), but indicating 
contribution from various redshift ranges. LL is adopted as the supernova model. 
These figures are taken from |17| . 



Figure Efb) shows the contribution by supernova neutrinos emitted from various 
redshift ranges. At high-energy region E u > 10 MeV, the dominant flux comes from 
the local supernovae (0 < z < 1), while the low-energy side is mainly contributed by 
the high- redshift events (z > 1). This is because the energy of neutrinos that were 
emitted from a supernova at redshift z is reduced by a factor of (1 + z)~ x reflecting 
the expansion of the universe, and therefore high-redshift supernovae only contribute 
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to low-energy flux. We also show the energy-integrated flux from each redshift range 
in table 121 in the case of the LL supernova model. From the table, it is found that in 
the energy range of our interest, more than 70% of the flux comes from local supernova 
explosions at z < 1, while the high-redshift (z > 2) supernova contribution is very small. 

3.2. Event rate at water Cerenkov detectors 

The water Cerenkov neutrino detectors have greatly succeeded in probing the properties 
of neutrinos as elementary particles, such as neutrino oscillation. The SK detector is one 
of these detectors, and its large fiducial volume (22.5 kton) might enable us to detect 
the diffuse background of SRNs. Furthermore, much larger water Cerenkov detectors 
such as HK and UNO are being planned. SRN detection is most likely with the inverse 
/3-decay reaction with protons in water, u e p — > e + n, and its cross section is precisely 
understood [HHl EH] • In our calculation, we use the trigger threshold of SK-I (before the 
accident). 

The expected event rates at such detectors are shown in figures 0a) and E^b) in 
units of /» (22.5 kton yr)" 1 MeV" 1 ; with SK, it takes a year to obtain the shown SRN 
spectrum, while with HK and UNO, much less time [1 yr x (22.5 kton/Vfid), where Vgd is 
the fiducial volume of HK or UNO] is necessary because of their larger fiducial volume. 
Figure Efa) compares the three models of the original supernova neutrino spectrum, 
and figure Efb) shows the contribution to the total event rate from each redshift range. 
In table El we summarize the event rate integrated over various energy ranges for three 
supernova models. The expected event rate is 0.97-2.3 /* (22.5 kton yr) -1 for E c > 10 
MeV and 0.25-1.0 /* (22.5 kton yr)- 1 for E e > 18 MeV. This clearly indicates that if the 
background events that hinder the detection are negligible, the SK has already reached 
the required sensitivity for detecting SRNs; with the future HK and UNO, a statistically 
significant discussion would be possible. This also shows that the current shortage of our 
knowledge concerning the original supernova neutrino spectrum and luminosity gives at 
least a factor of 2 [E v > 10 MeV) to 4 {E v > 18 MeV) uncertainty to the event rate 
at the high-energy range (actual detection range). We also summarize the contribution 
from each redshift range in the same table, especially for the calculation with the LL 
model. The bulk of the detected events will come from the local universe (z < 1), but 
the considerable flux is potentially attributed to the range 1 < z < 2. 

4. Detectability and observational upper limit 

4-1. Background events against detection 

In the previous section, we calculated the expected SRN spectrum at the water Cerenkov 
detectors on the Earth, but the actual detection is quite restricted because of the 
presence of other background events. In this paper, we follow a detailed consideration 
of these backgrounds by Ando et al |12j . There are atmospheric and solar neutrinos, 
antineutrinos from nuclear reactors, spallation products induced by cosmic-ray muons 
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Figure 3. (a) Event rate at water Cerenkov detectors in units of /» (22.5 kton yr) _1 
MeV -1 for three supernova models, (b) The same as (a), but indicating contribution 
from various redshift ranges. LL is adopted as the supernova model. These figures are 
taken from [T71 . 



Table 3. Event rate of supernova relic neutrinos. 



Event rate [/* (22.5 kton yr) _1 ] 



Model 


Redshift range 


E c > 10 MeV 


E c > 18 MeV 


LL 


Total 


2.3 


1.0 




< z < l a 


1.7 (77.5) 


0.9 (87.5) 




1< z < 2 a 


0.5 (20.6) 


0.1 (11.9) 




2 < z < 3 a 


0.0 ( 1.7) 


0.0 ( 0.5) 




3 < z < 4 a 


0.0 ( 0.1) 


0.0 ( 0.0) 




4 < z < 5 a 


0.0 ( 0.0) 


0.0 ( 0.0) 


TBP 


Total 


0.97 


0.25 


KRJ 


Total 


1.7 


0.53 



a Contribution from each redshift range to the total (0 < z < 5) value are shown in 
parentheses as percentages. 



and decay products of invisible muons. We show in figures EJa) and Efb) the flux and 
event rate of SRNs and these background events. 

The flux of the atmospheric neutrinos is usually calculated using MC method 
including various relevant effects (flux of primary cosmic rays, solar modulation, 
geomagnetic field, interaction of cosmic rays in the air, and so on), and in that simulation 
one-dimensional approximation is used, i.e., after the interaction of primary cosmic ray 
particles with air nuclei, all the particles are assumed to be moving along the line of 
the momentum vector of the primary cosmic ray particles [70]. There are many authors 
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Neutrino Energy [MeV] Positron Energy [MeV] 

Figure 4. (a) Number flux of SRNs compared with other background neutrinos, (b) 
Event rate of SRNs and invisible muon decay products. In both panels, LL is adopted 
as the original neutrino spectrum. These figures are taken from 12 . 



who calculated the atmospheric neutrino flux (see [H]-[ZS] for recent one- dimensional 
results). We use in this paper the flux calculated by Gaisser et al [66, 67J. More recently, 
results of three-dimensional flux calculations have been reported by several groups (see 
|74[ I7E] and references therein) . These calculations show slight increase of the flux at 
low-energy regions, although the flux below 100 MeV is not given. Therefore, it should 
be noted that there is possibility that the flux of atmospheric neutrinos in the relevant 
energy regime is higher than that adopted in this paper by about 20-30%. 

Solar neutrino flux is dominant at energy range below 19 MeV. We use the flux 
predicted by the standard solar model (SSM) in figure |U |7B] . Since the solar neutrinos 
are not v e but z/ e , the cross section for them is about two orders of magnitude smaller 
than that for z/ e . Furthermore recoil electrons scattered by solar neutrinos strongly 
concentrate to the opposite direction of the Sun, in contrast to the isotropic distribution 
of v e events. Therefore, the solar neutrinos are an avoidable background, not as critical 
as other events. At the same energy range corresponding to E u < 19 MeV, there is 
another serious background that becomes an obstacle also for solar neutrino detection, 
i.e., spallation products induced by cosmic ray muons. The event rate of the spallation 
background is several hundred per day per 22.5 kton, and it is extremely difficult to 
reject all of these events. Therefore the solar neutrino range [E v < 19 MeV), cannot be 
used as an energy window. 

The third background which we must consider is anti-neutrinos from nuclear 
reactors. In each nuclear reactor, almost all the power comes from the fissions of the 
four isotopes, 235 U (~ 75%), 238 U (~ 7%), 239 Pu (~ 15%) and 241 Pu (~ 3%) [77j. 
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Each isotope produces a unique electron anti-neutrino spectrum through the decay of 
its fission fragments and their daughters. The u e spectrum from 235 U, 239 Pu and 241 Pu 
can be derived using the semi-empirical formula with which we fit data of detected (3- 
spectrum from fission by thermal neutrons [THUZHl- ( 238 U undergoes only fast neutron 
fission and hence electron spectrum from 238 U cannot be measured by this kind of 
experiment.) Above 7 MeV, the number of (3 counts drops dramatically and fitting 
error becomes large. In addition, with this method, as we determine the maximum f3 
energy and derives the energy distribution below that energy, it is difficult to estimate 
the errors at the high-energy range. While the v e spectra in [THJ EHj are given as tables, 
we use for simplicity somewhat less accurate analytical approximation given in [80J. As 
a normalization factor we use energy- integrated u e flux at Kamioka, 1.34 x 10 6 cm -2 
s _1 , which are the summation of the flux from various nuclear reactors in Japan and 
Korea [ZZj - 

With these backgrounds we discussed above and from figure Ufa), we expect the 
energy window of SRN events ranging 19-30 MeV. However, electrons or positrons from 
invisible muons are the largest background in the energy window from 19 to 60 MeV. 
This invisible muon event is illustrated as follows. Atmospheric neutrinos produce muons 
by interaction with the nucleons (both free and bound) in the fiducial volume. If these 
muons are produced with energies below Cerenkov radiation threshold (kinetic energy 
less than 53 MeV), then they will not be detected ("invisible muons"), but successively 
produced electrons and positrons from the muon decay will be visible. Since this muon 
decay signal will mimic the z/ e p — > e + n process in SK, it is difficult to distinguish SRN 
from these events. The spectrum of this invisible muon events is shown in figure Efb), 
compared with the SRN spectrum. Therefore, even at the remaining candidate of the 
energy window 19-30 MeV, there is a huge background due to the invisible muons, 
resulting in no energy window for the SRN detection at present. 

4-2. Detectability at pure-water Cerenkov detectors 

We expect that SRNs can be most likely detected at water Cerenkov detectors such as 
SK, because the largest class of volume can be realized, enabling the most statistically 
significant discussions. The most serious problem is that there is no energy window as 
we described in the previous subsection. In the energy range 19-30 MeV, the SRN event 
number is estimated to be N$kn = 0.73/*(Kff/22.5 kton • yr) for the LL model, where 
we define the effective volume of detectors V e g by (fiducial volume) x (time) x (efficiency) , 
which is the relevant quantity representing the detector performance. On the other 
hand, the background level would be as large as N hg ~ 3.4(^/22.5 kton • yr), mainly 
contributed by the invisible muon decay products. Therefore, signal-to-noise ratio S/N 
can be written as 



S/N 
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plotted in figure as a function of /* for several values of effective volume. For detectors 
with the size of SK, one-year observation is found to be insufficient to obtain significant 
SRN detection. Further data-taking for about 10 years may reach S/N < a few, but 
this number is still insufficient in statistical significance. Future mega-ton class detectors 
such as HK and UNO now planned potentially give us a considerable number of SRN 
detection, by data-taking for a couple of years, as can be seen from the upper curve in 
figure El 




0.5 1 1.5 2 2.5 3 3.5 



Figure 5. Signal-to-noise ratio S/N of SRNs at pure-water Cerenkov detectors JSJ, 
as a function of a correction factor for the SFR model (0J. LL is assumed for the 
original neutrino spectrum. Each line is labeled by the value of the effective volume 
V cS . 



4-3. Detectability at other detectors 

In this subsection, we discuss the SRN detectability at SNO (the heavy water Cerenkov 
detector) and KamLAND (the scintillation detector), which is mainly discussed in 
[TTJ [12J and OHHj, respectively. An advantage of these detectors is that we are able to 
identify z/ e events using delayed coincidence signals; neutrons produced via z/ e d — > e + nn 
(SNO) and u e p — > e + n (KamLAND) are tagged by deuterons, CI or He (SNO) and by 
protons (KamLAND), resulting in 7-ray cascades, identified with the preceded positron 
signal. Using this technique, we can remove other backgrounds from non-z/ e origin (solar 
neutrinos, invisible muons and spallation products), opening up an energy window in 
the range of 10-30 MeV. Unfortunately, because the detector volume is small (both 
detectors are about 1 kton), the expected SRN rate is quite small, 0.03/* yr _1 for SNO 
and 0.1/* yr _1 for KamLAND [12] • More recently, Strigari et al ^B] estimated the 
event rate at KamLAND to be ~ 0.4 yr _1 , but they used the energy range above 6 
MeV, simply because there have been no events seen above this energy at KamLAND 
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[D3] . In addition, the detectability at liquid argon detectors has been discussed in [T^|; 
although the detection is still challenging, potential advantage is that these detectors 
can mainly capture u e , which is difficult with other detectors. 

4-4- Current observational limit and future prospects 

The most stringent upper limit on the SRN flux is obtained by the observation for 1496 
days (4.1 yr) at the SK detectors [20J. This limit is obtained by the statistical analysis 
including the background events from atmospheric neutrinos and invisible muons, and is 
< 1.2 cm~ 2 s _1 over the energy region of E v > 19.3 MeV (90% CL). Comparing with the 
prediction for the same energy range given in tabled we find that the current SK limit is 
only about a factor 2.5-8.5 larger than our prediction using the reference model for the 
cosmic SFR with /* = 1, depending on the adopted original neutrino spectrum. This 
strong constraint motivated many theoretical studies ^H] JD] an d has been translated 
into constraints on various quantities as we further discuss in the following sections. 

Although the SK limit was derived by a careful analysis of the spectral data over 
the energy range E e > 18 MeV, we here show that a very rough statistical argument 
using equation (JSJ) can fairly well reproduce the same limit. In equation (jSJ), we consider 
setting a limit at 90% CL, which corresponds to 1.64<r level (or S/N < 1.64), as a result 
of no SRN detection for V c s = 22.5 kton • 4.1 yr observation. Substituting these values 
into equation (jHJ), it is solved for /*, resulting in the upper limit /* < 2.6 (90% CL). 
This value is then translated into the flux limit using the relation between the SRN flux 
and /* (for the LL model) given in tabled i.e, F v < 1.2 cm~ 2 s _1 for E v > 19.3 MeV 
(90% CL), which is in remarkable agreement with the actual result reported by the SK 
group [2D] . 

In the near future, sensitivity of water Cerenkov detectors for the SRN detection 
would be significantly improved by the promising technique proposed recently [21]. The 
basic idea is the same as the delayed coincidence technique actually adopted by SNO 
or KamLAND (see discussions in section IP]) . but GdCl3 is dissolved into the pure- 
water of SK (or other future detectors), which enables us to actively identify u e by 
capturing neutrons produced by the u e p — > e + n reaction. Owing to this proposal, the 
range 10-30 MeV would be an energy window because we can positively distinguish 
the v e signal from other backgrounds such as solar neutrinos (z/ e ), invisible muon events 
and spallation products. The neutron capture efficiency by Gd is estimated to be 90% 
with the proposed 0.2% admixture by mass of GdCl 3 in water, and subsequently 8 MeV 
7-cascade occurs from the excited Gd. The single-electron energy equivalent to this 
cascade was found to be 3-8 MeV by a careful simulation [8*T] . and with the trigger 
threshold adopted in SK-I, only about 50% of such cascades can be detected actually. 
However, it is expected that SK-III, which will begin operation in mid-2006, will trigger 
at 100% efficiency above 3 MeV, with good trigger efficiency down to 2.5 MeV |21j . 
In that case most of the 7-cascades from Gd will be detected with their preceding 
signal of positrons. From this point on, we assume 100% efficiency; even if we abandon 
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Table 4. Constraints on supernova rate and star formation rate models. Upper part 
shows the SK limit on the SRN observation [201, whereas lower part shows the results 
of local supernova surveys and observational inference of the local SFR [53] • 





/* 


#sn(0) (yr- 1 Mpc" 3 ) 


V>*(0) (M yr- 1 Mpc- 3 ) 


LL 


< 2.6 


< 2.2 x 10~ 4 


< 1.8 x 10~ 2 


TBP 


< 8.6 


< 7.3 x 10~ 4 


< 6.0 x 10~ 2 


KRJ 


< 4.3 


< 3.6 x 10~ 4 


< 3.0 x 10~ 2 


□3 
M 


0.9-1.9 
0.7-4.2 


(0.8-1.6) xlO- 4 


(0.5-2. 9)xl0~ 2 



this assumption, it does not affect our physical conclusion, since the relevant quantity 
representing the detector performance is the effective volume V e s that already includes 
efficiency 

5. Implication for cosmic star formation history 

5.1. Constraints from the current observational limit 

The current SK upper limit on the SRN flux is already stringent to give some physical 
or astronomical consequences. In this subsection, we discuss a constraint on the cosmic 
SFR from the SRN limit; similar arguments have been given in jTHj. Comparing with our 
flux predictions using several supernova models (table |2J), the SK limit can be directly 
translated into the bound on the correction parameter /* to the cosmic SFR introduced 
in (j3J), since the redshift dependence is roughly consistent among various observations. 
Further, the constraint on /* can then be used as that on the local supernova rate 
and SFR. The results are summarized in table H] for various supernova models. For 
comparison, in the lower part of the same table, we also show the result of local supernova 
surveys, .Rsn(O) = (0.8-1.6) x 10 -4 yr -1 Mpc~ 3 [56J, and the observationally inferred 
SFR V*(0) = (0.5-2.9) x 10" 2 M yr" 1 Mpc" 3 0. In particular for the local SFR, 
the SK limit on the SRN flux may rule out some fraction of the observationally inferred 
value, if we choose the LL model for the original neutrino spectrum; for the other 
two supernova models, the SK limit is very close to the current upper bound on the 
cosmic SFR by observations with the light. Therefore, we stress that the neutrino 
observation has already reached sensitivity to the cosmic SFR comparable with the 
usual and traditional approaches using the light. 

5.2. Performance of Gd-loaded detectors 

Performance of the proposed Gd-loaded detectors as an SFR probe is of our interest, 
and has recently been investigated in detail by Ando ^2] using the MC simulation; in 
this subsection, we briefly introduce his discussion. Although we focus here on how 
far the SFR can be probed by SRN observation, the uncertainty from the supernova 
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neutrino spectrum would give a fair amount of error. However, this problem can be 
solved if a supernova explosion occurs in our galaxy; the expected event number is 
about 5000-10,000 at SK, when supernova neutrino burst occurs at 10 kpc, and it 
will enable a statistically significant discussion concerning the neutrino spectrum from 
supernova explosions. Even if there are no galactic supernovae in the near future, 
remarkable development of the supernova simulation can be expected with the growth 
of computational resources and numerical technique. With such developments, the 
supernova neutrino spectrum and luminosity may be uncovered, and the ambiguity is 
expected to be reduced significantly. Thus, in this paper we assume that the supernova 
neutrino spectrum is well understood and that our reference models are fairly good 
representatives of nature; we analyze the SFR alone with several free parameters. 

The basic procedure adopted in fH] is as follows. (1) The expected signal 
(spectrum) at a Gd-loaded detector in the range 10-30 MeV is simulated, assuming 
that there are no background events. In that process, the SFR given by equation (j^J) 
with /* = 1 and the LL model as neutrino spectrum are used for the generation of the 
SRN signal. (2) Then the SRN spectrum is analyzed using the maximum likelihood 
method with two free parameters of the SFR and a set of the best-fit values for those 
parameters is obtained; they are concerned with the supernova rate as 



where _Rg N represents the local supernova rate and a determines the slope of supernova 
rate evolution. Although the constant SFR at z > 1 is assumed, it is found that the 
result would be the same even if this assumption is changed. This is because the bulk 
of the detected event comes from local supernova as already shown. (3) 10 3 such MC 
simulations are performed and 10 3 independent sets of best-fit parameters are obtained. 
Then we discuss the standard deviation of the distributions of such best-fit parameter 
sets and the implications for the cosmic SFR. 

First we discuss the performance of Gd-SK for 5 years, or an effective volume of 
22.5 kton • 5 yr. Because the expected event number is only ~ 10, the parameters _Rg N 
and a cannot both be well determined at once. Therefore, one of those parameters 
should be fixed in advance at some value inferred from other observations. First, the 
value of i?g N was fixed to be 1.2 x 10~ 4 yr -1 Mpc -1 , which was inferred from the local 
supernova survey [HH|, and the distribution of the best-fit values of parameter a was 
obtained. The result of 10 3 MC simulations are shown in figure |Hf a) as a histogram of 
the distribution of best-fit parameters a (solid histogram) . The average value of these 
10 3 values for a is found to be 2.7, and the standard deviation is 0.8, i.e., a = 2.7 ± 0.8. 
Then in turn, the value of a was fixed to be 2.9 in order to obtain the distribution of 
best-fit values for the local supernova rate i?g N from the SRN observation. The result of 
10 3 MC generations and analyses in this case is shown in figure Efb). The average value 
for i?° N is 1.2 x 10~ 4 yr" 1 Mpc~ 3 , and the standard deviation is 0.4 x 10~ 4 yr -1 Mpc~ 3 , 
i.e., i?° N = (1.2 ± 0.4) x 10~ 4 yr -1 Mpc~ 3 . Figure [7| shows the comparison between the 




for z < 1, 
for z > 1, 



(9) 
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supernova rate model, in which the parameter is inferred from the MC simulations, and 
the "true" reference model; the cases of fixed i?g N and a are shown in figures [7(a) and 
UJb), respectively. The allowed region at the la level is located between the two dotted 
curves, while the solid curve represents the reference model, from which MC data were 
generated. Thus, with the Gd-SK detector we can roughly reproduce the supernova rate 
profile at z < 1 for 5 years operation. 




Figure 6. Distribution of 10 3 best-fit values for (a) a and (b) i?g N , which are obtained 
from the analyses of each MC generation. The effective volume is 22.5 kton • 5 yr 
for solid histogram and 440 kton • 5 yr for dashed histogram, (a) The value of the 
local supernova rate is fixed to be i?g N = 1.2 x 10~ 4 yr -1 Mpc -3 . The resulting 
distributions are characterized by a — 2.7±0.8 for an effective volume of 22.5 kton-5 yr 
and a = 2.5 ± 0.2 for 440 kton • 5 yr. (b) The value of a is fixed to be 2.9. The 
resulting distributions are characterized by i?g N = (1.2 ± 0.4) x 10~ 4 yr -1 Mpc~ 3 for 
22.5 kton • 5 yr and = (1.0 ± 0.1) x 10" 4 yr" 1 Mpc~ 3 for 440 kton • 5 yr. These 
figures are taken from |17|. 



Now we turn our attention to future mega-ton class detectors such as Gd-HK or 
Gd-UNO. With these detectors, the effective volume considered, 440 kton ■ 5 yr, is 
expected to be realized in several years from the start of their operation. First the 
same analysis as illustrated above was performed; i.e., one of relevant parameters, a 
or -Rg N , was fixed and the dependence on the remaining parameter was investigated. 
The result of these cases are also shown in figures Efa) and ijb) as dashed histograms, 
which give a = 2.5 ± 0.2 and i?° N = (1.0 ± 0.1) x 10" 4 yr -1 Mpc~ 3 , respectively. The 
statistical errors are considerably reduced compared with the case of 22.5 kton ■ 5 yr, 
because of the ~ 20 times larger effective volume. Thus, future maga-ton detectors 
will possibly pin down, within 10% statistical error, either the index of supernova rate 
evolution a or the local supernova rate i?g N if the other is known in advance. The 
dashed curves in figure UJ set the allowed region of the supernova rate at the la level 
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0.4 0.6 

Redshift z 

Figure 7. Supernova rate as a function of redshift. In both panels, solid curves 
represent the reference model, from which MC data were generated, (a) The allowed 
region at the la level, concerning the fitting parameter a with fixed i?g N , is shown as 
the area between the two dotted curves for an effective volume of 22.5 kton • 5 yr and 
as the area between the two dashed curves for an effective volume of 440 kton • 5 yr. 
(b) The same as (a) but for fitting parameter i?g N with fixed a. This figure is taken 
from Hlj. 



by the considered detectors, well reproducing the assumed model. In principle, we can 
determine both parameters by SRN observation, because i?g N is concerned with the 
absolute value of the flux alone but a is concerned with both the absolute value and the 
spectral shape; i.e., these two parameters are not degenerate with each other. Thus, the 
same procedure was repeated but without fixing the values of a or -Rg N . The distribution 
of 10 3 best-fit parameter sets of (a, -Rsn) was calculated assuming a detector with an 
effective volume of 440 kton • 5 yr; the mean values and the standard deviations were 
found to be a = 3.5 ± 1.3 and i? SN = (8.8 ± 4.8) x 10~ 5 yr" 1 Mpc~ 3 . Even though the 
effective volume is as large as 440 kton • 5 yr, it is still insufficient for determining both 
parameters at once. 

Throughout the above arguments or those in the reference SFR model 

with /* = 1 was assumed and they seem to be dependent on the adopted values 
of /*. However, the results can easily be applied to the case with other values of 
/*. This is because the statistical errors are inversely proportinal to the square root 
of the event number, which is simply proportional to f*V e g in our arguments; i.e., 
5a/(a), 5i2g N /(i?g N ) oc (/*Kfj) -1 / 2 . In consequence, we find that the discussions given 
above and in J7| are quite general and can be easily modified in the other cases of /* 
using the following relations: 

f V 1/2 5R ™~_ Q2 ( M* Y' 2 no) 

(a) ' V22.5 kton • 5 yr/ ' <i2g N > ' V 22 -5 kton • 5 yr y v ; 
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As a next step, it would be useful to investigate the possibility that how far (to 
which z) we can probe the cosmic SFR by the SRN observations. As we have already 
shown, the main contribution to the SRN even rate at 10-30 MeV comes from low- 
redshift region < z < 1. Signals from further high-redshift universe would become 
enhanced if we could reduce the lower energy threshold E t h of the SRN detection range. 
The value of _E t h is restricted to 10 MeV because at energy regions lower than this, 
there is a large background of reactor neutrinos; its removal is impossible with the 
current detection methods. Since the SK and HK detectors are (will be) located at 
Kamioka in Japan, they are seriously affected by background neutrinos from many 
nuclear reactors. If some large- volume detectors were built at a location free from such 
background, the lower threshold energy could be reduced, enabling us to probe the 
high-redshift supernova rate. Ando \Fj\ also discussed performance of future mega-ton 
class detectors as a function of the value of E t ^ and found that the behaviour of SFR 
at 1 < z < 2 more and more affects the detected event number at E th < E e < 30 MeV, 
as we reduce the threshold energy. 

6. Implication for neutrino properties 

6.1. Inverted mass hierarchy 

Throughout the above discussions, we have assumed normal hierarchy of neutrino 
masses (mi <C rriz). However, the case of inverted mass hierarchy has not been 
experimentally excluded yet, and we explore this possibility in this subsection following 
the discussion given by Ando and Sato j!4j . In this case, flavour conversions inside 
the supernova envelope change dramatically, compared with the normal mass hierarchy 
already discussed in section 12.41 Since z/3 is the lightest, z/ e are created as z/3, owing 
to large matter potential. In that case, it is well known that at a so-called resonance 
point, there occurs a level crossing between v\ and z/3 (for a more detailed discussion, 
see, e.g., |H2])- At this resonance point, complete v\ «-> z/3 conversion occurs when the 
so-called adiabaticity parameter is sufficiently small compared to unity (it is said that 
resonance is "non-adiabatic"), while conversion never occurs when it is large (adiabatic 
resonance). The adiabaticity parameter 7 is quite sensitive to the value of #13, i.e., 
7 oc sin 2 26 l i3; when sin 2 2^13 > 10 -3 (sin 2 2^13 < 10~ 5 ), the resonance is known to 
be completely adiabatic (non-adiabatic) [82J. When the resonance is completely non- 
adiabatic (because of small #13), the situation is the same as in the case of normal mass 
hierarchy already discussed in section 12.41 (because v c at production become v\ at the 
stellar surface), and the v e spectrum after oscillation is represented by equation (|7J). On 
the other hand, adiabatic resonance (due to large 6*13) forces v e at production to become 
z/ 3 when they escape from the stellar surface, and therefore the observed z/ e spectrum is 
given by 



dE, 
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Table 5. Flux of supernova relic neutrinos in the case of inverted mass hierarchy. 
Values given in this table are applicable only when the value of #13 is large enough to 
induce completely adiabatic resonance, i.e., sin 2 2^13 > 10~ 3 . If sin 2 26>i3 < 10~ 5 , on 
the other hand, the results become the same as those given in tabled 



Model 




Flux (/* cm 2 


s- 1 ) 


Total 


E v > 11.3 MeV 


E v > 19.3 MeV 


LL 


9.4 


3.1 


0.94 


TBP 


13.8 


1.9 


0.30 


KRJ 


12.7 


2.2 


0.38 



Table 6. Event rate of supernova relic neutrinos in the case of inverted mass hierarchy. 
Values given in this table are applicable only when the value of #13 is large enough to 
induce completely adiabatic resonance, i.e., sin 2 2^13 > 10 -3 . If sin 2 26*13 < 10~ 5 , on 
the other hand, the results become the same as those in table 01 



Model 


Event rate [/» (22.5 kton yr)" 1 ] 


E c > 10 MeV E c > 18 MeV 


LL 


3.8 2.3 


TBP 


1.6 0.58 


KRJ 


2.0 0.76 



The second equality follows from the fact that the value of | C/ e 3 1 2 is constrained to be 
much smaller than unity from reactor experiments jMj- Thus, equation (fTTj) indicates 
that complete conversion takes place between u e and u x . When the value of #13 is large 
enough to induce adiabatic resonance (sin 2#i3 > 1(T 3 ), the obtained SRN flux and 
spectrum should be very different from ones obtained in sections 13.11 and 13.21 The 
SRN flux and event rate in this case were calculated with equations © and (fTTjl . 
and the results are summarized in tables 03 and |H1 The values (with the LL model) 
shown in this table are consistent with the calculation by |TJJ, in which numerically 
calculated conversion probabilities were adopted with some specific oscillation parameter 
sets (which include a model with inverted mass hierarchy and sin 2 26 , 13 = 0.04), as well 
as realistic stellar density profiles; they also included the shock wave propagation and 
the Earth matter effect in their calculations, but both were found to affect only by a 
few percent and irrelevant. 

The total flux becomes 9.4-14 /* cm -2 s _1 , somewhat smaller than the values 
given m table because the total flux is dominated by the low-energy region. The 
fluxes at E v > 19.3 MeV are enhanced to be 0.30-0.94 /* cm -2 s _1 , but this is still 
below the current 90% CL upper limit of 1.2 cm -2 s -1 obtained by the SK observation 
if /* = 1 is adopted. The event rate at the future detectable energy range, E v > 10 
MeV, is expected to become 1.6-3.8 /*(22.5 kton yr) _1 , which is considerably larger 
than the values in the case of normal mass hierarchy, 0.97-2.3 /*(22.5 kton yr) _1 . The 
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increase (decrease) of the flux and event rate integrated over the high (total) energy 
region is due to the very high efficiency of the flavour conversion, z/ x — > u e , inside the 
supernova envelope; because the original is x are expected to be produced with larger 
average energy as shown in table ^ the efficient conversion makes the SRN spectrum 
harder, which enhances the flux and event rate at the high-energy region. Thus, if the 
inverted mass hierarchy, as well as the large value for # 13 , were realized in nature, SRN 
detection would be rather easier, compared with the other cases. Although we do not 
repeat the MC simulations that were introduced in section ElZJ the results can be easily 
inferred; the statistical errors in this case would be ~ (3.8/2.3) 1 / 2 = 1.3 times smaller 
than the values given in equation (fTU|) . because they are inversely proportional to the 
square root of the event number. 

6.2. Resonant spin- flavour conversion 

If neutrinos have non-zero magnetic moment as large as 10~ 12 /iB, where /ie is the Bohr 
magneton, it potentially changes supernova neutrino signal owing to an additional effect 
of the flavour conversions. Especially if neutrinos are the Majorana particles, the 
interaction between the Majorana magnetic moment and supernova magnetic fields 
induces a spin-flavour conversion (e.g., v e <-> v^ T ,v c <-> i/ M)T ) resonantly (see [EE] 
and references therein). This mechanism could potentially give quite a characteristic 
supernova neutrino signal at detectors, and also is expected to affect the SRN spectrum 
significantly. However, it is still premature to estimate the SRN flux including the spin- 
flavour conversions for several reasons. First, the shock wave propagation can change 
the magnetic field structure as well as the density profile of supernova progenitors, both 
of which are essential in calculating flavour- conversion probabilities [83]. Because there 
is no reliable supernova simulation that succeeded in pushing the shock wave outside 
the core, we are even not at the stage to start the calculation. Furthermore, calculating 
how the propagating shock wave changes the magnetic field structure would be quite a 
difficult task. The second reason is that the spin-flavour conversion probabilities strongly 
depend on the metallicity of progenitor stars [SI]. Because SRN is the accumulation 
of neutrinos from past supernova explosions, it should be quite natural that the poor 
metal stars give some contribution to the SRN flux. In consequence, it would be difficult 
to obtain some implication for the neutrino magnetic moment or supernova magnetic 
fields from the SRN detection at present; instead, a future galactic supernova neutrino 
burst might give some clues by its time profile or spectrum j^Sl EH] • 

6.3. Decaying neutrinos 

Non-radiative neutrino decay, which is not satisfactorily constrained, potentially and 
significantly changes the SRN flux and spectrum. Most stringent lower limit on the 
neutrino lifetime-to-mass ratio comes from the solar neutrino observations [SH] [II] as 
well as the meson decay experiments |I2] _ [II|, which is r/m > 10~ 4 s/eV. Since this 
limit is still very weak, neutrino decay could also change the detected signal from high- 
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energy astrophysical objects [13] or the galactic supernova explosions [HI]-[IE], and could 
alter usual discussions on the early universe and structure formation as well as on 
supernova coolings |S01^|S21- In this subsection, we discuss the potential consequence of 
the neutrino decay on the SRN flux, with which we can obtain the most stringent limit 
on the neutrino lifetime in principle, as pointed out by Ando 

Since the flux estimation of SRNs including the neutrino decay is somewhat 
complicated, we rather follow a simple and intuitive argument given in exact 
formulation was recently derived by Fogli et al [THj. Ando assumed the following 
conditions: (i) daughter neutrinos are active species; (ii) neutrino mass spectrum is 
quasi-degenerate (mi ~ vri2 ~ rn^); and neutrino mass hierarchy is normal. This is 
because the SRN flux would be most strongly enhanced owing to the decay and it might 
give a larger flux than the current observational limit by SK j2Hj; in that case the flux 
limit can be translated into the limit on the neutrino lifetime. He also restricted his 
discussion only to the helicity-conserving mode of z/3 — > i>\ and v<i — > V\ for simplicity, 
characterized by lifetimes r 3 and T2 of each mode; in models discussed in [HOJ — |HHJ , the 
above condition (ii) strongly suppresses the helicity-flip mode. In the actual calculation, 
instead of these lifetimes, he introduced the "decay redshift" zf (i = 2, 3) of the mass 
eigenstate i>i as two free parameters. If the source redshift z is larger than the decay 
redshift zf, all the neutrinos z/j were assumed to decay, on the other hand if z < zf, 
completely survive. The relation between zf and T%/m can be written as 



where we assume E v = 10 MeV as a typical neutrino energy when evaluating the lifetime 
Tj from zf. 

Figure |H1 shows the SRN flux for various parameter sets of decay redshifts (z$,zf) 
as a function of neutrino energy; the LL model as the original neutrino spectrum and 
f* — l were assumed. The solid curve in figure |H1 shows the SRN flux in the case that 
z 2 = z t = - 2 max(= 5), which represents the same flux as that without the neutrino 
decay. Then, the decay of the heaviest mass eigenstate z/ 3 was included by his reducing 
the value of zf, with keeping i> 2 stable. When zf = 1.0 (dotted curve), the SRN flux at 
low energy region (E v < 35 MeV) deviates from the no-decay model. This is because the 
neutrinos from supernovae at redshift larger than zf = 1.0 are affected by the z/3 — > v\ 
decay and it results in the increase of i> e . Since the neutrino energies are redshifted 
by a factor of (1 + z)" 1 owing to an expansion of the universe, the decay effect can 
be seen at low energy alone. When the value of zf is reduced to 10~ 2 , the neutrinos 
even from the nearby sources are influenced by the z/ 3 — > V\ decay, resulting in the 
deviation over the entire energy range as shown by the long-dashed curve in figure El If 
the v-i — > v\ decay is added, it further enhances the SRN flux. In table [TJ the SRN flux 
integrated over the energy range of E v > 19.3 MeV is summarized, for the each decay 
model. In the second column the lifetime-to-mass ratio is indicated, which corresponds 
to each decay redshift, which is obtained using equation (JT2J). The ratio between the 
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Table 7. Predicted SRN flux for various decay models. The LL model is assumed as 
the original neutrino spectrum and integrated energy range is E v > 19.3 MeV. The 
ratio of the prediction and the limit is shown in the fourth column. 





(T 2 /m,T 3 /m) 


Predicted flux 




Model (4.4) 


(s/cV) 


(/» cm" 2 s- 1 ) 


Prediction/limit 


(5.0,5.0) 


(3.9 x 10 10 ,3.9 x 10 10 ) 


0.43 


0.35/, 


(5.0,1-0) 


(3.9 x 10 10 ,2.4 x 10 10 ) 


0.55 


0.42/* 


(5.0,0.2) 


(3.9 x 10 10 ,7.7 x 10 9 ) 


0.93 


0.75/* 


(5.0,10- 2 ) 


(3.9 x 10 10 ,4.4 x 10 8 ) 


1.0 


0.88/, 


(10-2,10-2) 


(4.4 x 10 8 ,4.4 x 10 8 ) 


1.4 


1-2/* 



predicted flux and the 90% CL upper limit given by the SK observation [20] is also 
shown in the fourth column. This result shows that several decaying models may have 
already been excluded or severely constrained by the current SRN limit by SK. It also 
proves that the SRN observation can potentially give us the lower limit on the neutrino 
lifetime as strong as ~ 10 10 s/eV, which is many orders of magnitude stronger than that 
by the solar neutrino observations. In addition, future detectors such as the HK and 
UNO (maybe loaded with Gd) is expected to greatly improve our knowledge of the SRN 
spectral shape as well as its flux; if a number of data were actually acquired by such 
detectors, the most general and model-independent discussions concerning the neutrino 
decay would become accessible. 




10 20 30 40 50 60 

Neutrino Energy [MeV] 

Figure 8. SRN flux for various parameter sets of decay redshifts, calculated with the 
LL model and /* = 1. Each label represents (z%, z$). This figure is taken from |15j . 
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7. Conclusions 

As stressed several times in this paper, we are at an exciting era of the SRN 
observation. This is because the current SK upper limit [20] is just above several 
theoretical predictions using realistic models for the cosmic SFR and supernova neutrino 
spectrum, and also because of the promising technique proposed in [21] could greatly 
improve performance of the presently working or planned future detectors, especially 
for detecting SRNs. Being at such a situation, we believe that it is timely to review 
recent progress of theoretical and observational researches concerning SRNs from various 
points of view. Many physical and astronomical consequences can be derived from the 
constraints on the SRNs as we have reviewed through this paper. 

Basics of SRN calculation were given in sections |2] and H3 Models involved in 
the SRN calculations are those of the cosmic SFR and original supernova neutrino 
spectrum. Although the SFR-z relation has been intensively studied in recent years, 
there remains a fair amount of uncertainty; in fact, the best estimate of the local SFR 
density ranges fairly widely as (0.5-2.9) xl0~ 2 hi yr _1 Mpc -3 54J. In order to take 
this uncertainty into account, we have introduced the correction factor /* in equation 
(jlj), for which a wide range (0.7-4.2) is still allowed from the SFR observation. As 
for the original supernova spectrum, contrary to the traditional approach using the 
Fermi-Dirac distribution, we followed the approach taken by our group [121 Ej, i.e., 
adopting the results of numerical simulations. Especially following [T7| . we used three 
neutrino spectra by different groups, i.e., LL [5Z], TBP [HE] and KRJ |09J, and performed 
comparison among models. Recent progress of neutrino experiments has proven that 
neutrinos are massive and mix among different flavours. This effect was also taken into 
account appropriately. As the result, it was found that the uncertainty concerning the 
original neutrino spectrum gives difference in the resulting SRN flux by at least a factor 
of 2-4. The expected event rate at a water Cerenkov detector with a fiducial volume 
of 22.5 kton was found to be 0.25-1 /* yr" 1 (for E c > 19.3 MeV) and 1-2 /* yr" 1 (for 
E c > 10 MeV), depending on these supernova models (see table EJ). 

Besides the flux estimation, a careful discussion about background events is 
definitely necessary in order to investigate the detectability. This has been thoroughly 
argued in ^2) and we followed them in section EJ The most serious background comes 
from decay products of invisible muons, which are produced almost at rest inside the 
detector and that is why "invisible." Because of the background, the detection is severely 
restricted; for a pure-water Cerenkov detector with the size of SK, it would take about 
(maybe more than) 10 years in order to reach the signal-to-noise ratio of a few (see 
figureEJ). We also discussed performance of other detectors such as SNO and KamLAND 
in the same section, by mainly following [TSl ll6|. The current observational upper limit 
by SK [2Hj and the proposed technique that potentially improves the performance of 
water Cerenkov detectors [ST] were briefly introduced there. In addition, we showed 
that the SK limit on the SRN flux can be fairly well reproduced by a simple statistical 
argument using equation (jHJ). 
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Section was devoted to the discussion on the current status and future prospects 
for obtaining limits on the cosmic SFR from the SRN observations. We first showed 
that the current SRN limit from SK (2H| already sets a stringent limit on the local SFR 
value, as discussed in JT3 . Surprisingly, it might rule out a part of local SFR value 
allowed by the current observations using the light as shown in table EJ Performance 
of future detectors loaded with Gd as proposed in [2T] was investigated in detail using 
the MC simulation, following discussion in ^7]. The distribution of the best-fit values 
of two parameters concerning SFR (jHJ) was obtained by the generated 10 3 MC data 
and by accompanying likelihood analyses. If one of two free parameters was fixed 
somehow in advance using other observations, then the SRN observation can constrain 
the remaining parameter by accuracy of ~ 0.3(/*Kff/22.5 kton ■ 5 yi)' l/2 \ this number 
would soon become quite significant if the future mega-ton class detectors such as HK 
or UNO loaded with Gd started data-taking. 

Finally in section 13 SRN constraints on particle physics models of neutrinos were 
discussed. Neutrino oscillations in the case of inverted mass hierarchy FJ3], resonant spin- 
flavour conversions induced by the neutrino magnetic moment (see jHH] and references 
therein) and neutrino decay CHj are expected to considerably change the SRN signal 
at the detectors. In particular for neutrino decay model, the SRN observation potentially 
gives the most stringent constraint on the neutrino lifetime, compared with a galactic 
supernova [IH]-|IH|, high-energy astrophysical objects [IS] as well as the most reliable 
(but weak) limit at present by the solar neutrino observations |39j-|41j. 
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